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We construct a model of an excitable medium with elastic rather than the usual diffusive coupling. We 
explore the dynamics of elastic excitable media, which we find to be dominated by low dimensional structures, 
including global oscillations, period-doubled pacemakers, and propagating fronts. We suggest that examples 
of elastic excitable media are to be found in such diverse physical systems as Burridge-Knopoff models of 
frictional sliding, electronic transmission lines, and active optical waveguides. 



PACS numbers: 91.30.-f, 83.50.By, 84.40.Az, 87.22.As 

The response to a perturbation is the defining characteristic 
of an excitable element. Above a certain threshold in ampli- 
tude, a perturbation will excite a quiescent element that then 
decays back to quiescence after a characteristic time insen- 
sitive to the magnitude of the perturbation, during which it 
is unresponsive to further perturbations. Excitable elements, 
when coupled to their neighbors into an assembly, become an 
excitable medium. Models of excitable media with diffusive 
coupling, such as those of van der Pol, FitzHugh, and Nagumo 
|Jl|] have successfully described pattern formation phenomena 
in biology and chemistry, and have also captured the atten- 
tion of many working in the area of nonlinear science because 
of their complex spatiotemporal dynamics [0]. On the other 
hand, the properties of excitable media with elastic rather than 
diffusive coupling have not to date been investigated. We 
present here some physical systems — frictional sliding, elec- 
tronic transmission lines, and active optical waveguides — as 
diverse examples of such elastic excitable media, and ana- 
lyze their dynamics, which we find to be dominated by low 
dimensional structures, including global oscillations, period- 
doubled pacemakers, and propagating fronts. We investigate 
the stability of the oscillations, and estimate the velocity of 
the fronts. 

Whilst touching on electronic and optical applications of 
elastic excitable media, we shall concentrate on frictional slid- 
ing. The Burridge-Knopoff model [[|] was originally intro- 
duced as a representation of earthquake fault dynamics. It 
describes the interaction of two tectonic plates in a geological 
fault as a chain of blocks elastically coupled together and to 
one of the plates, and subject to a friction force by the surface 
of the other plate, such that they perform stick-slip motion. In 
this Letter, we show that a Burridge-Knopoff model in which 
the usual Carlson-Langer velocity weakening dry-friction law 
|Q] is replaced by the original Burridge-Knopoff lubricated 
creep-slip version showing viscous properties at both the low 
and high velocity limits, is an elastic excitable medium. 

The reintroduction of this friction law in the Burridge- 
Knopoff model is motivated by experiments showing that it 
represents friction in a range of materials, including paper on 
paper metal on metal [Q], and rock on rock M\, and by the- 
oretical arguments M\. The model well represents the qualita- 



tive characteristics of laboratory stick-slip experiments with 
elastic gels [0], Further interest derives from studies of peel- 
ing adhesive tape jl0[], of Saffman-Taylor fracture in viscous 
fingering and of the Portevin-Le Chatelier effect of dis- 
continuous yielding Jl2|], showing that stick-slip phenomena 
in those systems are induced by this same form of friction law. 

Burridge & Knopoff introduced a class of simple models 
that describe the contact region between two tectonic plates 
as a chain of N blocks of equal mass to, mutually coupled by 
springs of Hooke constant k c and equilibrium length a. The 
blocks are pulled by the bulk of one plate moving at velocity 
V through constant elastic shear k p against the friction Ff 
between the two plates, as shown in Fig. [j](a). In the stationary 
frame, the equation of motion for the ith block is 

mXi = k c (X l+1 - 2X, + AVi) 

-kpiXi-VQ-FfiXi), (1) 

where Xi is the departure of block i from its equilibrium po- 
sition. Usually, following Carlson & Langer [Q], the friction 
is taken to be asymptotically velocity weakening, so that the 
blocks perform stick-slip motion: any individual block sticks 
to the surface until the pulling force exceeds the static friction 
threshold. Once the block starts slipping, the dynamic fric- 
tion diminishes monotonically with its velocity, as Fig. |l](b) 
illustrates. While this multivalued form of the friction law is 
perfectly admissible in the discrete model of Eq. (jl|), it clearly 
poses a problem if we attempt to take the continuum limit 




FIG. 1. (a) The Burridge-Knopoff model. (b) The ve- 
locity weakening stick-slip friction law of Carlson & Langer 
[Ff(v) — Fq sgn(ii)/(l + \v\), where v is the velocity of the block] 
(solid line) and the Burridge-Knopoff type creep-slip friction law 
we use (dashed line), showing the threshold (vertical line) where the 
block starts to slip. 



1 



To overcome this difficulty, one may replace the static fric- 
tion discontinuity by a small region of high viscosity. Phys- 
ically, one can see this cut-off as representing lubrication ef- 
fects that produce stable creep at low velocities, making our 
model creep-slip rather than stick-slip. Since we wish to 
model the frictional characteristics displayed by a range of 
materials [§-0], we also have the friction become viscous (ve- 
locity strengthening) at high velocities, rather than decay to 
zero as in the Carlson-Langer version. Our friction model 
(Fig. [jjb)), which is similar to that originally introduced by 
Burridge & Knopoff, thus comprises three regions: first ve- 
locity strengthening, then velocity weakening, and finally ve- 
locity strengthening again. We shall refer to this model as 
asymptotically velocity strengthening, to be contrasted with 
the asymptotically velocity weakening friction usually consid- 
ered. With this in mind, the continuum limit of Eq. (|l|) in 
dimensionless variables (see, e.g., [pj|]) is 

X = c 2 x"-(x-^)-7<Mx)- (2) 

x{x,t) represents the time-dependent local longitudinal de- 
formation of the surface of the upper plate in the static refer- 
ence frame of the lower plate; dots and primes are temporal 
and spatial derivatives, respectively. In this continuum limit, 
the number of blocks N becomes the system size S. 4>(x) 
is our normalized asymptotically velocity strengthening fric- 
tion. There are three dimensionless parameters: 7 measures 
the magnitude of the friction, c is the longitudinal speed of 
sound, and v represents the pulling velocity or slip rate. From 
Eq. (||), we obtain an expression for the local velocity ip = x 
of the interface that, written as a couple of differential equa- 
tions of first order in time, gives us our continuum Burridge- 
Knopoff model 

ip = 7(V-<t>W), (3a) 
f] = -- f - 1 {ip-v-c 2 ip"). (3b) 

With the proper choice of the function (j)(ip), and if we 
may disregard for a moment the positioning of the spatial 
derivatives term, Eqs. (Q) constitute a version of the van 
der Pol-FitzHugh-Nagumo model that is the prototypical de- 
scription of an excitable medium. 

While this type of medium was first studied in the physi- 
ology of cardiac and nervous tissues, the same properties are 
also seen in chemical and physical systems. An electrical car- 
icature of an element of an excitable biological membrane is 
given by the circuit shown in Fig. ||(a). The membrane rep- 
resented by the capacitor C is charged by ion pumps char- 
acterized by the current generator i and drained through a 
nonlinear resistance ni across the membrane. This nonlinear 
element should have the v — i characteristic tj> shown in Fig. 

, and could be implemented by a tunnel diode, a neon 
lamp, or any electronic device with a range of negative resis- 
tance. The inductance L models the finite switching time of 
the ion channels in the membrane. The circuit equations are 
then formally equivalent to Eqs. (^J) above, where ip and r\ are 
proportional to the current and potential difference across the 



nonlinear device, respectively, and the current io is propor- 
tional to v. Such circuit elements have been used in electronic 
experiments with Burridge-Knopoff models [OJ. In the mod- 
els of excitable media of FitzHugh & Nagumo, <p is generally 
taken to be as in Fig. ^(b): it thus has exactly the same form 
as used in our Burridge-Knopoff model (cf. Fig. |j](b)). 
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FIG. 2. An electronic excitable medium, (a) An element of the 
medium: the circuit and (b) the v — i characteristic of ni. (c) Resistive 
(diffusive) spatial coupling, (d) Capacitive (elastic) spatial coupling 
— an active transmission line. 

The spatial derivatives in Eqs. (||) are unusual in the realm 
of excitable media. In a biological excitable membrane, the 
excitation propagates in space by diffusion. With the electri- 
cal model described above one can build a spatially extended 
diffusive medium by coupling resistively several excitable el- 
ements, as indicated in Fig. ^(c). Such a coupling would give 
rise to current diffusion: the term involving ip" would appear 
in Eq. ( |3a| ) instead of in Eq. (|3b"|). As they are written how- 
ever, Eqs. (||) represent the network of capacitively coupled 
excitable elements illustrated in Fig. ^(d). This is a discrete 
representation of an active transmission line, where L and C 
are the distributed inductance and capacitance and ni is the 
distributed gain achieved by a nonlinear negative resistance 
along the line. The distributed serial capacitance C c repre- 
sents the effects of dispersion in the line. The continuum limit 
of this network is a caricature of an active optical waveguide 
such as a fiber amplifier or laser. 

Further motivation for the study of Eqs. (j^) is provided by 
laboratory experiments [^| in which an elastic gel was placed 
in a Taylor-Couette type apparatus. Stick-slip events oc- 
curred at the contact surface of the gel with the rotating inner 
cylinder. In a large part of the parameter space, low dimen- 
sional structures — continuous slip, global oscillations, and 
propagating fronts — were found rather than complex spa- 
tiotemporal patterns that might have been expected to dom- 
inate. With this experiment in mind we focus on periodic 
boundary conditions. Since we are interested in qualitative 
aspects of the dynamics, we concentrate on the case where <p 
is taken to be a cubic polynomial. By shifting the origin of 
the variables ip and 77 and the parameter v we can reduce <p to 
the form 4>(ip) — ip 3 /3 — ip. This makes our model a set of 
elastically coupled van der Pol oscillators. For this election of 
<f>, the slip threshold, defined as the crossover as we increase v 
from velocity strengthening creep to velocity weakening slip, 
lies at v = —1. 
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FIG. 3. Spatiotemporal plots of ij){x,t) where time is vertical, 
space is horizontal, (a) A plot of the solution with stable uniform 
global oscillations for v = —0.92. (b) Stable propagating fronts at 
v = —0.96. In both cases, the other parameter values are 7 = 2 and 
c = 0.3. The system size S = 20, and the time elapsed 50. 

The continuous slip solution, where one surface moves uni- 
formly with respect to the other, was examined by Carlson & 
Langer [pj[]. They showed that the growth rate of perturba- 
tions of any wavelength at a particular slip rate is minus the 
slope of the friction function at that slip rate. This implies 
that with asymptotically velocity weakening friction the con- 
tinuous slip solution is always unstable, while in our case this 
solution is stable for high and low slip rates where the slope is 
positive. Global oscillations, where one surface moves period- 
ically in time with respect to the other, are also stable at some 
slip rates in our model. Figure ||(a) shows the spatiotempo- 
ral pattern generated from an arbitrary initial condition with 
the system relatively far above the slip threshold. Here all 
points of space oscillate in phase with an instantaneous ve- 
locity that follows the dynamics of the van der Pol oscillator. 
To confirm the stability of these oscillations, we have calcu- 
lated their Floquet multipliers. Global oscillations are solu- 
tions of Eqs. ([|) with the spatial derivative set to zero, which 
correspond to the period-T limit cycle of the van der Pol os- 
cillator **(t + T) = $*(t) where = (ip*(t),7]*(t)), 
which is stable against spatially homogeneous perturbations 
for — 1 < v < 1 (i.e., above the slip threshold). We consider 
small spatial perturbations of the form \&* (t) + e[t) exp(iqx). 
Dropping quadratic and higher powers of e(t) in Eqs. (g), 
the resulting linear equations for e(t) have time-periodic co- 
efficients, and by the Floquet theorem their general solution 
is e{t) = £oP(t) exp(A q t), where P(t) is a period-T ma- 
trix and A q a time-independent matrix. Once the limit cy- 
cle has been numerically determined, the matrix exp(A g T) 
may be computed by numerically integrating the linearized 
equations over the period T. Its two eigenvalues, the Floquet 
multipliers (roj , m q 2 ^), which may be either real, or com- 
plex conjugates, determine the stability of global oscillations 
under perturbations of wavenumber q. The growth rate per 
period of such a perturbation is the maximal Floquet expo- 
nent X q = maxi{K[ln(mq*' ) )]}. The dispersion relation X q 
as a function of q relatively far above the slip threshold is 
never positive, implying stability under perturbations of all 
scales. This contrasts with the same calculations performed 



for asymptotically velocity weakening friction, when there is 
always a range of long wavelengths for which perturbations 
grow exponentially, so that global oscillations are always un- 
stable when this wavelength is within the system size. This 
is a significant difference between our model and Burridge- 
Knopoff models with asymptotically velocity weakening fric- 
tion. Such global oscillations have been noted in laboratory 
friction experiments with elastic gels [0]. 





FIG. 4. (a) A family of dispersion relations for the Floquet expo- 
nent showing the destabilization close to the slip threshold of global 
oscillations against perturbations of a single wavenumber q c ~ 3. 
Shown are dispersion relations for (from the lower to the upper 
curve) v = -0.90, -0.92, and -0.94. (b) Spatiotemporal plot of 
i]){x, t) for v = —0.94 showing the development of period-doubled 
structures resulting from this destabilization. Time elapsed is 125, 
other settings are as in Fig. |5| The spatial periodicity is determined 
by Qc- 

For slip rates not far above the slip threshold, global os- 
cillations become unstable. The dispersion relation in Fig. 
|](a) shows that the destabilization occurs at a single finite 
wavenumber q c . Furthermore, the Floquet multiplier m 9c as- 
sociated with the maximal Floquet exponent crosses the com- 
plex unit circle through —1 at this point, indicating a period- 
doubling bifurcation. Figure |](b) is a spatiotemporal plot 
showing the nature of the structure arising near this insta- 
bility. Small perturbations from the global uniformly oscil- 
lating regime grow to become structures with a well-defined 
wavenumber q c . For system size S, such structures evolve into 
a configuration of q c S /(2tt) synchronized pacemakers emit- 
ting fronts in both directions. These fronts annihilate each 
other at the corresponding q c S/ (2ir) points, that then in turn 
become pacemakers. The resulting spatiotemporal structure 
is such that at an arbitrary point in space the periodicity in 
time is twice that of the original oscillations. These period- 
doubled structures are not well known in pattern formation, 
so their presence here is an interesting theoretical prediction. 

At slip rates right above the slip threshold, the front solu- 
tions cease to annihilate each other and evolve into fronts that 
propagate right around the system (Fig. [|(b)). One can get an 
analytical handle on these propagating fronts by positing so- 
lutions of the type -0(a;, t) = f(z), where z = x/v + t, and 
v is the front velocity. This ansatz with the further rescaling 
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7 1 — c 2 /v 2 leads to 



4 + / 



(4) 



which is again the van der Pol equation, but with the non- 
linearity rescaled by fj, — 7/yl — c 2 /v 2 . The propagating 
fronts are then periodic solutions of the van der Pol equa- 
tion. The parameter fi is undefined until the value of the 
front velocity v is chosen. However, we know that the pe- 
riod of the solution is a function T = T(fi) of /x: in the 
limit of large /1, T behaves as T = kfj, + Of/i -1 ), where 
k = 3 + {v 2 - l)ln[(4 - v 2 )/{l - v 2 )] [[jj. Since this 
period should be commensurate with the system size S, we 
have the condition nT(fi(v)) — S/ (vy/l — c 2 /v 2 ), where n 
is an integer, to select the allowed front velocities, which in the 
large [i limit gives us the quantizing condition v = S/ (nk-f). 
Because Eq. (|]) has bounded solutions only if v 2 > c 2 , the 
propagating fronts are supersonic. Such supersonic propagat- 
ing fronts have also been noted in Burridge-Knopoff models 
with asymptotically velocity weakening friction [[j~6j[l7|]. 

The original Burridge-Knopoff model was introduced as a 
means to reproduce the gross features of the statistics of real 
earthquakes [|[|. The Gutenberg-Richter power-law distribu- 
tions [ |l8| | in the statistics of slip events in the model have been 
considered as an example of self-organized criticality and 
have been related to the presence of infinitely many degrees 
of freedom in the system fll^]. However, recent numerical ex- 
periments [|l6 20] indicate that power-law distributions of slip 
events may be due to discretization, finite size, and transient 
effects, and are not present at long times in the continuum 
limit. This has lead to a questioning of the relevance of the 
model to earthquakes in the real world. Our aim here, how- 
ever, has been to investigate a Burridge-Knopoff model with 
a physically interesting friction law relevant to laboratory fric- 
tion experiments: ours are simply laboratory earthquakes. 

We have returned to the asymptotically velocity strength- 
ening type of friction law originally introduced by Burridge 
& Knopoff, which was abandoned by Carlson & Langer and 
later investigators, and with it have been successful in ex- 
plaining some results of laboratory friction experiments. Our 
model constitutes a form of excitable medium, not previously 
considered, with elastic instead of diffusive coupling between 
spatial elements. We have studied the spatiotemporal dynam- 
ics of elastic excitable media and find novel dynamics. We 
believe that this type of elestic excitable medium may have 



applications beyond laboratory friction experiments to elec- 
tronic transmission lines and active optical waveguides. 
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